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We introduce the formalism of generalized Fourier transforms in the context of risk management. 
We develop a general framework to efficiently compute the most popular risk measures, Value- 
at-Risk and Expected Shortfall (also known as Conditional Value-at-Risk). The only ingredient 
required by our approach is the knowledge of the characteristic function describing the financial 
data in use. This allows to extend risk analysis to those non-Gaussian models defined in the Fourier 
space, such as Levy noise driven processes and stochastic volatility models. We test our analytical 
results on data sets coming from various financial indexes, finding that our predictions outperform 
those provided by the standard Log-Normal dynamics and are in remarkable agreement with those 
of the benchmark historical approach. 
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INTRODUCTION 



September 2008 financial crisis has dramatically high- 
lighted the need for reliable, easy to understand and im- 
plement instruments to measure and manage risk. The 
high volatility of financial markets during the Nineties 
induced academics and practitioners to design sophisti- 
cated risk management tools. According to the recently 
revised capital adequacy framework, commonly known 
as Basel II accord [l|, any financial institution has to 
meet stringent capital requirements in order to cover the 
various sources of risk to be faced as a result of normal 
operations. Today the most widely used measure to man- 
age market risk in the financial industry is Value-at-Risk 
(VaR). VaR refers to the maximum potential loss over 
a given period at a certain confidence level (CL) and 
can be used to measure the risk of individual assets and 
portfolios of assets as well. Because of its conceptual sim- 
plicity, VaR has been extensively adopted by regulators 
and it generally provides a reasonably accurate estimate 
of risk. However, VaR is known to suffer from important 
drawbacks: it can violate the sub-additivity rule for port- 
folio risk, which is a required property for any consistent 
measure of risk 0, Q , and it does not quantify the typical 
loss incurred when the risk threshold is exceeded. The 
expected shortfall (ES), defined as the expected loss con- 
ditional on the VaR threshold being exceeded, overcomes 
these disadvantages and leads to more consistent results. 

Three main approaches are known in the literature and 
used in practice to compute VaR and ES: the parametric 
approach, the historical one, and Monte Carlo simula- 
tions of the stochastic dynamics of a given stock price re- 
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turns model. The parametric approach usually relies on 
the normality assumption for the returns distribution, al- 
though some analytical results using non-Gaussian func- 
tional forms are available in the literature [1, I3j S @]- 
However, it is well known that empirical price returns ex- 
hibit heavy tails and a certain degree of asymmetry; the 
historical simulation approach is often used in order to 
capture their leptokurtic nature. The last approach con- 
sists of Monte Carlo simulations of the return dynamics, 
but it usually requires very intensive simulations to get 
to acceptably accurate risk estimates. As a result of the 
present situation, reliable and hopefully fast methods to 
calculate financial risk are mandatory. 

In this article we shall present a general framework 
to compute VaR and ES by only relying on the knowl- 
edge of the closed-form characteristic function (CF) de- 
scribing the distribution of the financial returns under 
analysis. Our approach draws on the original ideas de- 
veloped in 0, [1] by A. L. Lewis and A. Lipton, who 
introduced generalized Fourier calculus in the context of 
option pricing, and extends them to the risk management 
framework. The advantage of our approach is manifold. 
The equations we obtain for risk measures are intuitive 
and easy to read, since both VaR and ES turn out to be 
expressed in terms of one same function evaluated for dif- 
ferent arguments. Moreover, the evaluation of such for- 
mulae is computationally efficient, since running twice a 
fast Fourier transform (EFT) algorithm yields both VaR 
and ES values over wide ranges of the CL. This is a re- 
markable feature, since under standard approaches risk 
measures would need to be recomputed every time the 
CL is changed. Fourier inversion based approaches can 
be found in literature- the first attempt dates back to 
the work of Rouvinez [9|, and later developments are dis- 
cussed in [Tol. [Til [T^. [l3l H. [Tsj . However, the Fourier in- 
version is usually employed to compute an approximation 



2 



of the cumulative function through the Inversion Theo- 
rem [l^ , then the quantile corresponding to the fixed CL 
is computed by root-finding algorithms. However, this fi- 
nal step has to be iterated over the entire set of desired 
CLs, while our approach directly provides the risk esti- 
mates. A very useful and easy to interpret graphical rep- 
resentation of the results can also be sketched. Finally, 
being based on the use of CFs, our method is readily ap- 
plicable to a number of interesting distributions whose 
probability density function (PDF) is not known analyt- 
ically. Remarkable examples are represented by the class 
of Levy distributions, both in their original and truncated 
versions. 

In this article we focus on Truncated Levy Distribu- 
tions (TLDs), which, given their ability to reproduce 
some of the stylized facts observed in real market data, 
have been introduced and applied in the context of fi- 
nancial analysis by several authors [3 Holl . The 
approach we propose can also deal with stochastic volatil- 
ity models (SVMs), which have already been successfully 
used in the context of derivative pricing. Interestingly, 
to the best of our knowledge, these models have never 
been employed in risk management before. Our frame- 
work is naturally suited to models that are well defined 
in terms of the CF such as the Stein-Stein [HI HI], Hes- 
ton [2^, [2^ , Schobel-ZhuJ^ , and exponential Ornstein- 
Uhlenbeck [2^, [H, [23] models. We choose to work 
with the Heston model, whose popul arity is rapidly grow- 
ing amongst financial practitioners [30|. A further fea- 
ture of SVMs we wish to investigate in the present work 
is their ability to provide high order normalized cumu- 
lants with different time scalings w.r.t. those implied by 
the Central Limit Theorem (CLT). We plan to test their 
performances when projecting risk estimates over time. 

The paper is organized as follows: in Section II the 
technical definitions of VaR and ES are recalled, and their 
expressions in terms of the CF are derived. In Section III 
the models we use to test our approach are presented. In 
Section IV the fitting procedures and the data analysis 
we performed are described, and the numerical results 
obtained for the risk measures are detailed. In Section V 
some conclusions are drawn. 



II. FORMULAE FOR RISK ESTIMATION 

The Value-at-Risk (VaR) is defined as the maximum 
potential loss —A* (with A* > 0) not to be exceeded 
at a given significance level P* £ (0, 1) over a fixed time 
horizon At. Thus, when considering the price variation 
AS = S — So oi some asset or index, VaR is implicitly 
defined by the integral equation 



-A* 



d(AS') p{AS) 



(1) 



fiAt, with fi the linear returns mean, equation ([T]) can be 
rewritten as 



-L* 



da; p(x) 



(2) 



where p is the PDF associated with x and L* = — ln(l — 
A*/S'o)-l-/iAt. It is worth mentioning that A* represents 
VaR in monetary units; we shall often use the normal- 
ized VaR A* = A* /So, which is usually presented as a 
percentage quantity (percentage VaR). 

In order to derive new expressions for VaR, we adapt 
the approach developed by A. L. Lewis and, indepen- 
dently, by A. Lipton in the context of derivative pricing 
under stochastic volatility 0, 0. Let us represent the 
PDF p in terms of the generalized Fourier transform / 



p{x) 



2tt 



+oo+ii/ 



deb f{4>) e-'*^ 



In this integral expression (/) = uj + iiy is a complex vari- 
able (w, v £ M), whose imaginary part v belongs to the 
proper strip of regularity (i/_, z/+) of the extended char- 
acteristic function (ECF) /. Such a strip is delimited 
by the possible singularities of / (which can be shown 
to be purely imaginary under suitable conditions 7]) ly- 
ing closest to the origin in the complex upper (lower) 
half plane, whose imaginary part reads (^^-)- With 
these positions, we can plug the previous expression in 
equation ^ and switch the integration order. Thus, we 
obtain 



P* = 



1 

2^ 



+oc+ii/ 



-L* 



dx e""^^ 



and to ensure convergence of the second integral we re- 
quire V to be strictly positive. So, restricting v e (0, 
equation ([2]) eventually becomes 



27r 



-\-oc-\-ii/ 



d0 /(</>)- 



2tt 



duj 



f{uj + iiy) e} 



iuL* 



-vL' 



-Re 



did ■ 



where the final equality is obtained by exploiting the sym- 
metries of the real and imaginary parts of the ECF. Then, 
defining the function 



r+00 

G^(L*,6i) = e-^^' / duj 
Jo 



allows us to write 



Re Gy{L*,v) 



-So 



(3) 



(4) 



where p is the PDF describing AS'. When switching to 
the centered logarithmic returns x = ln(l -I- AS/ So) — 



We now follow a similar line of reasoning for the Ex- 
pected Shortfall (ES), defined as the average potential 
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loss when the VaR threshold for a fixed P* is exceeded. 
With the same notation as above, in terms of the linear 
and centered logarithmic returns, we can write for the 
ES E* 



E*{V*) = 

_ 5*0 



_A*(P*) 
-So 

-L*(P*) 



d(A5) AS p{AS) 



dx p{x) (e^+'^^* - 1) , (5) 



where we have made explicit the dependence of E* , A* 
and L* on P*. From now on we shall drop this depen- 
dence and we shall assume 5*0 = 1. As we did before, we 
can plug the generalized Fourier transform of p into the 
previous equation and switch the integration order 



E* = — 



-\-oo-\-\u 

( 



27rP' 

X / dx (e^i-'-^)^ 



J 



(6) 



When we pose (f) = uj + iv, the first of the two integrands 
in dec requires v > —\ to he evaluated, while the sec- 
ond one requires v > Q. So, all in all, we are again left 
with V £ (0,1^+) and, recalling the definition ([3]) of G^, 
equation ([6]) reads 



E* ^1 



^^A. ReG.(£*,z.+ l) 



(7) 



Formulae ([4]) and ([7|) do represent the first main con- 
tribution of this work and let us now see what the main 
advantage in their use is. Usually, VaR and ES are evalu- 
ated in correspondence of a single fixed significance level 
P* by means of different techniques, see [2| and the on- 
line repository at www.gloriEmiundi.org for an exhaus- 
tive review. In the best case scenario, the financial prac- 



titioner is provided with a closed-form expression depen- 
dent on distributional assumptions and parametric in the 
quantile of a standardized PDF and in few free parame- 
ters to be calibrated on the financial time series in use, 
see e.g. [1, [Sll, [S^. However, every time the value of 
P* is changed, risk measures need to be re-computed. 
Moreover, as anticipated in the introduction, up to now 
there are no efficient ways to conjugate risk estimation 
with models fully characterized in terms of the CF. In- 
deed, this is the case for the class of Levy stable distribu- 
tions and their exponentially damped version [itI . fisj and 
for all those dynamical models emerging in the context 
of option pricing under stochastic volatility, such as the 
Stein-Stein, Heston, Schobel-Zhu, exponential Ornstein- 
Uhlenbeck models iH, [H, HI, IH, [l?, '28^ and their 
extensions dealing with jump diffusion, e.g. [33] . The use 
of the G,y function can overcome all these drawbacks. As 
a matter of fact, once the ECF of the financial dynamics 
at hand is known in closed-form, a grid of to values can 
be set and {L* , v) (for an admissible value of i^) can 
be efficiently evaluated via FFT algorithms. This leaves 
us with a vector L* , which can be easily converted into a 
vector A* of VaR estimates. So, inserting L* in equation 
(|4|, we are quickly and efficiently provided with the full 
relation between the VaR estimates and the correspond- 
ing appropriate significance levels. Analogously, the FFT 
computation of G,y{L*,iy + 1) provides us with the E* 
spectrum over a whole range of significance levels. Thus, 
equations ^ and ^ provide a global information about 
the VaR and ES distributions over a wide range of P* 
values. These results lead to a very intuitive and easy to 
read graphical representation that we shall discuss in the 
final section. 



Let us now mention that the real part of the G^, func- 
tion can be put into a different form. In fact, we can 
define the function 



J 



/,(L*,0) =Re G,(L*,( 



-t-oo 



6*2 



;J^cos{ojL*)[9 Re/(, 



LU + iv) — UJ 



Imfiuj + iiy)]} 



-t-oo 



duj 



J I sin(wi*) [w Re/(w + 'w) + 6 Im/(w + ij/)] | (8) 
I 



and we can consequently rewrite the equations for risk 
measures as 



P* = 



E* = l-e 



(9) 



Now, the Ii, function in ([8]), being the sum of sine and 
cosine transforms, is perfectly suited to numerical evalua- 
tion by means of trapezoidal integration algorithms, and 
this partially prevents the risk estimates in ^ from be- 



ing affected by the FFT approximations, which, in turn, 
can be made negligible only by setting a large u) grid. 
The equations in ^ may prove to be very helpful to 
the evaluation of risk measures at a few specific values 
of the significance level P* with a very high precision. It 
is also worth stressing again that, remarkably, the only 
input those equations require is the ECF of the financial 
dynamics under study. As we shall discuss in the next 
section, this fact makes it possible to introduce and sue- 
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cessfully employ a number of models in the framework of 
risk management. 



III. MODELS 

One of the crucial points in risk analysis is the evalu- 
ation of risk measures over different time lags. So, gen- 
erally, a projection over the desired time horizons of the 
PDFs employed to model financial returns has to be per- 
formed. The most interesting case, being the one re- 
quired by regulators [H, is to project from one to ten 
trading days. Equations ([U, ([7]) and Q provide a very 
natural framework where different time scaling behaviors 
can be compared. In this section we discuss those aris- 
ing when considering two class of models: the first one 
of purely additive processes governed by the CLT and a 
second class of SVMs. 

As it has already been stressed, the main ingredient 
required by our approach is the knowledge of the closed- 
form ECF associated with the model in use. As a rep- 
resentative member of the class of additive process, we 
consider a very simple one given by an arithmetic mo- 
tion whose driving noise is described by a TLD XI} ■ The 
most appealing feature of the TLD is the ability to re- 
produce some of the stylized facts commonly observed in 
financial markets, such as the asymmetry and the excess 
of kurtosis [13]. As it is well known, the time scaling 
is governed by the CLT and this causes high order nor- 



malized cumulants to decrease monotonically with time. 
The scaling behavior is governed by a power law whose 
characteristic exponent is specific for the order of the cu- 
mulant, e.g. the skewness scales as while the kur- 

tosis scales according to the law. On the other hand, 
SVMs are naturally provided with a different, exponen- 
tially damped time scalings, not necessarily monotonic, 
which could be able to better capture real market data 
time scalings. In this article we want to test the ability 
of these models to capture the projection over horizon 
of risk measures in comparison with historical estimates 
and the standard Log-Normal dynamics. Since the most 
commonly used SVMs are well characterized in terms of 
the CF, formulae ([4]) and ^ allow for a proper compar- 
ison and the extension of the risk analysis to the context 
of these models does represent the second main contribu- 
tion of our work. In the following we focus on the Heston 
model [23] , since it represents a benchmark model in the 
option pricing framework. However, our app roach could 
also be easily applied to other SVMs [H [gH, [H, and 
extended models dealing, for example, with jump diffu- 
sion 13311. 



A. Truncated Levy Distributions 

The CF of a TLD can be expressed as f{(j>) = 
exp{H{(f))), where the cumulant generating function, or 
Hamiltonian, H is given by [H, [s^ 



J 



m) = - 



A2-7 



2 7(1-7) 



{(l+/?)exp[7log(A + i(/))] +(l-/3)exp [7 log(A - i(/.)] - 2A^} 



(10) 



and this gives rise to the following asymptotic behavior 
for the PDF ptl- 



Ptl{x) 



\x\ 



++00 e" 



-\\x\ 



[1-f /3sign(x)]. (11) 



,7,A,/3 is a constant depending on the four free param- 
eters which define the distribution. As it is clear from 
PT|) . (3 determines the level of asymmetry of the PDF, 
7 e (0, 2] is the tail exponent (7 = 2 reducing to the Nor- 
mal case), while A > is the decay factor; E > defines 
the level of the second moment of the distribution. The 
values of such parameters single out a particular TLD 
univocally. However, in our case, a corrective positional 
term has to be added to the Hamiltonian function PH)) . 
This is because we work with centered log-returns empir- 
ical distribution and so we need to correctly center the 
TLD on the real data. Thus we consider the modified 
Hamiltonian H ((/>) 



(12) 



where ki = — i 



dH 

d0 



0=0 



The singularity of the model rel- 



evant for the strip of regularity is readily found by solving 
the equation uj + iv+ = iA. 

The presence of an exponential damping in (jlip guaran- 
tees the finiteness of the variance, which can be shown 
to be equal to S^, and of all the higher order moments. 
The CLT applies and this has immediate consequences 
on the time scaling of many relevant quantities. As a 
matter of fact, upon the addition of N independent iden- 
tically distributed {i.i.d.) TLD variables the generic n-th 
order cumulant fc„ scales linearly with N. Then, for the 
skewness C and kurtosis k we have 



CiN) 



'1/2 



k{N) 



kl{N) 



N- 



(13) 

As it will be discussed in the next paragraph, the Heston 
stochastic volatility model leads to a much richer time 
evolution of the cumulants. 
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B. Heston dynamics 

The Heston model is defined by the two following cou- 
pled stochastic differential equations (SDEs) 

dS{t) = iiS{t)dt+ ^/^S{t)dWi{t) (14) 
dv{t) = a{a^ -v(t))dt + ky/v(t)dW2it) (15) 

under the initial conditions S{Q) = So and v{Q) — a^, 
and the Feller condition 2aa^ > k^. The noise W2{t) 
is defined in terms of the standard Brownian increments 



dWi and dW{t) through the usual relation dW2{t) = 
pdWi{t) + a/1 - p^dW{t) and p e [-1; 1]. In the follow- 
ing, we shall work with centered log-returns X, whose 
SDE can be derived from the previous ones: 

dX{t) = -^dt + ,/^dWiit) (16) 

with Xq — 0. By means of standard techniques [l^, [2^ . 
|36[ the cumulant generating function of the model can 
be shown to be 



aa 



{m - vW)t - 2 In (1 - g((/.)e-''(*)*) + 2 In (l - .g(0)) 



<^'^{^W~V{4>)) l-e-^W* 
1 - 5(0)e-''('^)* 



(17) 



with 



Data Sets and Calibration 



i]{(f>) 



a — ipk(j} 



From PT)) the cumulants of the Heston model PDF can 
be derived explicitly, and their analytical expressions 
are reported in Appendix |^ The identification of the 
strip (0,i^+) is more tricky for the Heston case. The 
relevant singular points solve the equations ri{ii''^) = 

0. Restricting to the case 



and 1 - .g(iz/^)e-''('''+)* 
P e ( 
2ap - 



1,1), from the former equation we have 
k+ ^yi2ap 



I [2fc(l 



ky + 4q;2(1 - p2 

the latter can not be solved explicitely but, once a set 
of parameters values has been fixed, we can numerically 



check if there is any positive Vj^ < 



However, for 



tipical values of a ~ 10 and k ~ 10, and for ly e (0, 1] 
we obtain positive t] and ^ of order 10^; so, for each 
t > 0, a candidate has to be strictly greater than one. 
Thus, if we fix j/ = 1 when integrating equations ^ and 
([7]), we have only to check that = i/"^ > 1. 



IV. DATA ANALYSIS AND NUMERICAL 
RESULTS 



As already remarked, our aim is to capture the empir- 
ical scaling of the returns and to exploit it in order to 
end up with reliable risk estimates projections over the 
time horizons of interest. Both in the TLD and in the 
Heston model case, we perform this task by focusing on 
the calibration on the first four cumulants time scalings. 



TLDs are fitted on empirical distributions by means of 
a simple step-by-step procedure according to which the 
free parameters in (jlOp . S, 7, /3 and A, are evaluated one 
at a time. First, extending the approaches developed in 
[itIIT^ to the case of asymmetric TLDs, the tail exponent 
7 is fitted exploiting the time scaling of the empirical zero 
return probability. Actually, the quantity pl(A' = 0) of a 
Levy distribution (a good approximation for the central 
region of a TLD, since the exponential damping mainly 
affects the tail regions) is used. Such a point can be 
shown to scale linearly with time t on a log-log scale 



logPL^(O) = log 



/(7,/3) 

CTT"/ 



-\ogt 

1 



(18) 



where c is a constant, while / is a function of the 7 and (3 
parameters. Clearly, from the slope of this linear relation 
the value of 7 can be estimated. Moreover, we consider 
the following relations [s^ valid for the variance, skew- 
ness and kurtosis of a TLD over an horizon t 



k2{t) = 

Kit) = (2-7)(3-7)/(A^S2i) 
at) = /3(2-7)/(aVe5^). 



(19) 



S can be readily estimated from the first relation, while 
A and /3 can be obtained from the remaining ones. Both 
last two relations can be rearranged into linear ones on 
a log - log scale. Thus, the parameters of interest can all 
be estimated via Marquardt-Levenberg linear fits. 

The Heston model calibration is performed differently. 
By imposing the initial condition vq — <t'^, we assumed 
the model to be in its stationary volatility state. The 
empirical mean p is estimated from the linear returns 
directly. Moreover, it is clear from the analytical expres- 
sion reported in Appendix [XI that a can be fitted on the 
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time scaling of the first cumulant. This leaves us with 
three more parameters to be estimated, i.e. a, p and k. 
We obtain the optimal values by solving the following 
minimization problem numerically 



10 



,k* = 



argmm 



I,l);a,fe>0^-^1 ,=2 



kfijAt)-kiijAt) 



efe.(jAt) 



In the previous equation the first sum runs over trad- 
ing days {At = 3.98 x 10~^ years), while the second one 
runs over the cumulants: kf represents the ?— th empiri- 
cally estimated cumulant, with statistical uncertainty ek^, 
whereas ki represents the i— th analytical cumulant. 

The calibration has been performed on three differ- 
ent data sets, made of 5000 daily returns each, from the 
German DAX 30 Index (from November 14th 1988 to 
September 9th 2008), the French CAC 40 Index (from 
November 10th 1988 to September 9th 2008) and the 
Dow Jones EURO STOXX 50 (SX5E) Index (from March 
10th 1989 to September 9th 2008). In Table [Jail of the 
relevant parameter estimates, both for the TLD and the 
Heston model, are reported. In Figure [1] the different 
time scalings obtained from such parameters are com- 
pared with the empirical ones. In the two upper figures 
the time scaling of the mean and variance is considered. 
An excellent agreement is observed for both the TLD 
and the Heston model. In the lower figures, even though 
the calibration was performed over the cumulants, we 
report the skewness and kurtosis time scaling because 
of their major relevance in risk analysis. Due to error 
propagation, the error bars look much more irregular in 
these cases. As it can be seen, the Heston model better 
describes the empirical scaling than the TLD model, es- 
pecially for the kurtosis data. For the skewness scaling, 
Heston slightly outperforms the TLD approach, as con- 
firmed in Figure [2] for the SX5E Index case. We shall see 
in the next section how this behaviour is translated in 
terms of the risk estimates. 



B. Risk Estimates 

In this section we detail and discuss all of the risk esti- 
mates we obtained. We present the results based on the 
TLD and Heston models under the generalized Fourier 
transform approach. As already discussed, this essen- 
tially amounts to the use of formulae ([3]), ([4]) and ([7]). 
As benchmark models we consider the Log-Normal dy- 
namics and the historical approach. The former is meant 
to provide a comparison with the results obtained under 
the standard normality assumption for the returns; the 
latter, instead, provides some insight into the actual risk 
levels of an asset. Historical estimates are obtained with 
standard methodologies (see [1] for example). 

We also employ the calibration procedure described in 
the above section in a bootstrap framework in order to 
provide the risk estimates with 68% CL intervals. For 
each data set we generate Mb = 1000 synthetic copies 



of the original time series and this is done by means of a 
GARCH(1,1) model to preserve the correlation structure 
of the volatility. Such a model is defined by the following 
couple of equations 



C- 
K 



crtzt 



Aaf 



(20) 



where Yt is the log-return at time t, at describes its 
volatility and the zt's, often referred to as innovations, 
correspond to a Gaussian white noise; C, K, G and 
A are constant quantities. The model calibration can 
be succesfuUy performed with the help of R software 
(www.r-project.org). Being z.z.d. Gaussian variables, 
the bootstrap technique can be applied to the innova- 
tions in order to generate replicas of the original time 
series preserving the volatility clustering. In our analy- 
sis, the TLD and Heston model calibrations have been 
carried out on each bootstrap copy. This provides us 
with a different set of parameter values for each copy 
that can be plugged into the ECF and consequently into 
equations (g]) and ([7]) to obtain copy-dependent risk esti- 
mates 9* for j = 1, . ■ . , Mb, where 9 can represent either 
VaR or ES. Then, bootstrap confidence levels are defined 
as [0* ; 9l_^], with the boundaries of the interval satisfy- 
ing the following relation: Prob(0* < 6** j_^) = a, 1 — a. 
Thus, a 68% CL interval requires a = 16%. The boot- 
strap technique allows us to draw statistically robust con- 
clusions. In Figure [3] we report two example bootstrap 
histograms of AJ^ and E'^ for the CAC 40 Index; from 
their analysis we can identify the extremes of the confi- 
dence intervals at the desired CL. 

As already mentioned, one of the main motivations to 
our approach is the need for a natural projection over 
horizon framework. As required by regulators, we fo- 
cus on risk measures evaluation at 1 day and 10 days 
horizons. Besides this, VaR and ES are evaluated at 
P* = 1%, this also being requested by regulators, and at 
P* = 5%. In Table HI] and Table the risk estimates 
for P* = 1% at 1 day and 10 days horizons are detailed; 
the "Hist" subscript refers to historical estimates, while 
N refers to the Normal ones, T and H to the TLD and 
Heston models, respectively. A few comments need now 
to be made. First of all, the best agreement with histor- 
ical estimates is found for AJ^ and E^, and it is remark- 
able for the P* = 1%, 1 day VaR. The agreement slightly 
worsens for the TLD, while the Normal estimates widely 
underestimate A^-^^ and i?Hist ■ These results confirm the 
ability of both the Heston and TLD models to better de- 
scribe tail events of the empirical distributions than the 
Log-Normal model. When considering a 10 days horizon 
the errors get much larger, since the empirical risk es- 
timates are evaluated on the basis of 500 returns only. 
The Heston and historical values are in best statistical 
agreement, while TLD and Normal estimates definitely 
worsen, and they get much closer as a consequence of 
the CLT, as expected. These results suggest that the 
projection over time horizon associated with the dynam- 
ics of the Heston model can provide a better description 
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TABLE I: Values of /j, (left) and of the TLD (center) and Heston (right) models parameters as estimated from the data sets. 
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FIG. 1: Comparison between the empirical, TLD and Heston time scaling of the mean, variance, skewness and kurtosis for the 
CAC 40 Index time series over a 10 days horizon. 
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FIG. 2: Comparison between the empirical, TLD and Heston time scaling of the skewness and kurtosis for the SX5E Index 
over a 10 days horizon. 
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TABLE II: Historical, Normal, TLD and Heston VaR and ES TABLE III: Historical, Normal, TLD and Heston VaR and 
estimates at P* — 1% over a 1 day horizon. ES estimates at P* = 1% over a 10 days horizon. 
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of the risk level for low P* than the one induced by the 
CLT. As far as a higher P* level is concerned, we consider 
P* = 5% (see Table IV] and Table |V|. One can see that, 
both for the daily and 10 days horizon, the difference be- 
tween the different VaR estimates reduces. For the ES, 
which by definition exibits a higher sensitivity to the tail 
behaviour than VaR, we again find the best performance 



for the Heston model. In particular, when switching to 
the 10 days horizon, the Heston model definitely leads 
to the best overall agreement with historical estimates. 
Finally, we present a graphical representation, very effec- 
tive for practical applications. In Figure [51 left panel, we 
plot A* against the significance level over a daily horizon; 
the curves are obtained via an adaptive trapezoidal inte- 
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TABLE IV: Historical, Normal, TLD and Heston VaR and 
ES estimates at P* = 5% over a 1 day horizon. 





AHist (%) 


(%) 


A*T (%) 


Al, (%) 


DAX 
CAC 
SX5E 


9 1 o+O.UU 


2 

^■^^-0.07 
1 nn+0.07 


2 01+"-"** 
1 oQ+o.oe 

'°-0.06 


9 90+0.13 
^•^O-O.IO 

2.08iH? 
2.011?-! 




Ekst (%) 


(%) 


E't (%) 


i^H (%) 


DAX 
CAC 
SX5E 


o nn+0.10 
'^■""-O.IO 

2 q4+0-i° 

^■^^-0.11 


9 yc+U.US 
^- ' '-'-0.09 
9 f,o+0.07 

9 1-1+0.09 
^■^^-0.09 


2.95±;!:- 

9 yi+0.20 
^- ' ^-0.19 


17+0.24 

■^•^ ' -0.24 


TABLE V: Historical, Normal, TLD and Heston VaR and ES 


estimates at P* = 5% 


over a 10 days 


horizon. 
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gration of the I,^{L* , 9) function to increase the numeri- 
cal accuracy w.r.t. FFT based approaches. We consider 
a grid of one hundred equaUy spaced P* values ranging 
from 0.1% to 10% and v = \. Triangles represent histor- 
ical estimates and the matching with the Heston curve is 
quite evident. We explicitly draw a vertical dotted line 
for P* = 1%. The crossing point with the Heston curve 
identifies the VaR estimate for the specified significance 
level. If we translate the estimate over the ES curve we 
obtain a new crossing point, whose projection over the 
horizontal axis returns the desired ES value. In the right 
panel of the same figure the analysis is performed for 
a 10 days horizon. This graphical approach shows how 
straightforward it is to obtain the risk estimates and to 
compare performances provided by different models. 

V. CONCLUSIONS 

In this paper we have shown how the extension of the 
CF to the complex domain can be successfully employed 
to derive compact formulae describing VaR and ES, the 
two market risk measures most financial industries ordi- 
narily use. The same technique has already been widely 



adopted and tested in the context of option pricing un- 
der stochastic volatility. The analogy is not so surpris- 
ing when noticing that the integral equation linking L* 
and P*, see equation ([5]), clearly resembles the relation 
between the strike and the price of an option with a 
digital payoff. A similar argument applies to the ES, 
even though in this case the payoff function looks a little 
bit more complicated, i.e. [1 - 9(a; -I- L*)] [e^+t^^*- - I), 
where Q is the Heaviside step function, see equation ([5J). 
Exploiting these analogies, we have obtained new integral 
representation of risk in terms of the function G^{L* , 6), 
particularly suited to efficient numerical integration using 
FFT algorithms. Based on the function Iu{L* , 6) and on 
adaptive trapezoidal algorithms, we have also suggested 
an alternative approach to perform an accurate integra- 
tion of our formulae. Our focus was on two types of log- 
returns stochastic dynamics. The first one is a simple 
arithmetic motion whose random increments correspond 
to a Truncated Levy noise, while the second is the Hes- 
ton stochastic volatility dynamics. Both models are an- 
alytically well-defined in terms of the CF; moreover, the 
CLT applies to the former predicting a power-law scaling 
with time of high order normalized cumulants. In the 
Heston case, instead, their time evolution obeys an ex- 
ponentially damped scaling. Since the risk measure pro- 
jection over time horizons is one of the points addressed 
by regulators, we have compared the performances pro- 
vided by these conceptually quite different models. We 
have tested our analytical formulae on a data set of fi- 
nancial indexes, the German DAX 30 Index, the French 
CAC 40 and the European Dow Jones EURO STOXX 
50. The results we obtained both for the TLD and He- 
ston cases have shown an excellent agreement with the 
historical benchmark values, within the bootstrap 68% 
CL . In particular, at the significance level required by 
regulators, i.e. P* = 1%, the Heston model was found 
to provide the best risk characterization, at least for the 
data we took into account. 

Possible perspectives of the present work concern the 
application of the approach to the evaluation of risk mea- 
sures using other SVMs, such as exponential Ornstcin- 
Uhlenbeck models and their extensions with jump diffu- 
sion. 



APPENDIX A: HESTON MODEL CUMULANTS 

In the following we report the analytical expressions of 
the Heston model cumulants 
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fcl = 

k2 = 
fcs = 

= 



— - \ - fc2e-2"* + Ake-'^Hk - 2ap) + 2at{4.a^ - k^ - Aakp) + k(8ap - 3k)] , 

^jfc^e"^"* - 3afce-2«* [-kt{k - 2ap) - 2{a - kp)] - Be""* [2akt{k - 2apf + 3fc^ - 8a^p - l&ak'^p 
+8a2fc(l + 2/92)] ^ ^_^3 _^ 8Q,3p ^ gQ,^2^ _ 4q,2^(-i ^ 2p2)] + gfc^ - 2Aa^p - A2ak'^p + &a^k{2, + 8p''')|, 
3fcV2 



{ - 3fc''e-4"* - 8fc2e-3"* [2QA;i(fc - 2ap) + 4q2 + fc2 _ g^^^j _ ^^-2at 40,2)^2^2 _ 2^,^) 



64q!7 

+2afct [A;^ - 16a^p - 12afcV + 4afc(3 + 4^^)] + 8a^ - 3fc^ - 32a^A;p + Safc^p + ISa^fc^p^ 
-Se""* - 2a^kt^{k - 2apf - 8at[k'^ - lak^p + Aa^p^ - 8a^kp{l + p^) + a'^k'^{3 + lAp^)] - 9k^ + IQak^p 
+S2a^kp{A + 3p2) - 16a''(l + Ap^) - Aa^k'^{9 + AOp^)] + 4at [5fc^ - 40aA;^p - 32a^fc/9(4 + 3/9^) + 16a''(l + Ap"^) 
+24aV^(l + V)] - 7^k^ + 5AAak^p + 12Sa^kp{7 + 6p^) - 32a^(3 + 16p^) - 6Aa^k^{A + 1 
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